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Abstract 
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Two conceptually different approaches (broadband-based ERBE and narrowband-based 
ISCCP approaches), used to derive the TOA longwave radiative fluxes, are compared using the 
ScaRaB simultaneous narrowband and broadband measurements. This study directly shows that 
the ERBE MLE-derived cloud covers implicitly contain some information on the cloud optical 
properties. A spurious view-zenith-angle dependence of the MLE scene identification scheme is 
confirmed by this study. Except for very thin cirrus clouds, differences between the ERBE and 
ISCCP approaches are in general < 10 W m' 2 for the TOA LW radiative fluxes. For clear pixels, 
the model calculated (ISCCP approach) TOA LW radiances are systematically smaller than the 
observations. Though the bias is found to be correlated on the column precipitable water amount, 
the exact source of this discrepancy remains undetermined and merits further study. Compared 
with the radiative transfer model used in this study, the ERBE LW ADMs are too weakly limb- 
darkened for optically thin clouds, but too strongly limb-darkened for optically thick clouds, 
indicating that more accurate instantaneous TOA LW flux estimations from the ERBE approach 
would require additional cloud classes based on cloud height and optical thickness. 
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1. Introduction 

The delicate balance between the incoming shortwave radiation and the outgoing longwave 
radiation at the top of the atmosphere (TOA) and its subtle regional variations is a direct indicator 
of the net gain/loss of radiant energy by our planet that determines the state of the atmospheric 
general circulation and the climate. Accurate knowledge of these two flux components is crucial 
for understanding the climate system and predicting climate change due to human activities, since 
the atmospheric circulation is forced by the vertical and horizontal gradients (i.e., differences) in 
radiative and latent heating/cooling, not by their magnitude. The TOA outgoing longwave 
radiation (OLR) can be directly estimated from satellite measurements of broadband radiances, 
such as the Nimbus-7 Earth Radiation Budget (Nimbus-7 ERB; Jacobowitz et al. 1984), the Earth 
Radiation Budget Experiment (ERBE; Barkstrom and Smith 1986; Barkstrom et al. 1989), the 
Scanner for Radiation Budget (ScaRaB; Kandel et al. 1994), and the ongoing Clouds and the 
Earth’s Radiant Energy System (CERES; Wielicki et al. 1996). Alternatively, OLR can also be 
indirectly inferred from narrowband radiances, either by a regression-based narrowband-to- 
broadband conversion technique (e.g., Ohring et al., 1984; Gruber and Krueger 1984) or by a 
detailed radiative transfer calculation with the cloud and surface radiative properties retrieved 
from narrowband radiances (e.g., the International Satellite Cloud Climatology Project, ISCCP; 
Schiffer and Rossow 1983; Rossow and Schiffer 1991) and other correlative atmospheric and 
surface datasets as input to the radiative transfer model (e.g., Rossow and Lacis 1990; Zhang et al. 
1995; Rossow and Zhang 1995; Hatzianastassiou and Vardavas 1999; Chen et al. 2000a, b). The 
narrowband-to-broadband conversion technique is still used in some operational earth radiation 
budget (ERB) products at present (e.g., Gruber and Krueger 1984) despite the limitations in the 
single-channel narrowband to broadband algorithm (Gruber et al. 1994). With the advent of more 
accurate global observations of cloud, atmosphere and surface properties from satellites, the 
approach which calculates broadband fluxes from retrieved physical quantities becomes more 
attractive since it intrinsically separates the effects of the clouds on the radiation field from other 
atmospheric and surface factors so that more direct estimates of the radiative effects of clouds can 
be made. Moreover, the modeling approach also allows us to partition the TOA radiation budget 
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into its surface and in-atmosphere components in a self-consistent manner (Rossow and Lacis 
1990). 

Comparison between the model calculated OLR from ISCCP retrieved cloud properties, 
supplemented by other atmospheric and surface datasets, and the ERBE OLR shows that while the 
two fluxes agree reasonably well in general, systematic differences exist in specific locations 
suggesting particular local surface, atmospheric and cloud conditions are responsible for the 
disagreement between these two products (Rossow and Zhang 1995; Xu 1997). It is imperative to 
identify the sources of these discrepancies, because these systematic differences, though relatively 
small in magnitude, are certainly significant when compared with the radiation budget changes 
that appear interannually or that are being considered in relation to possible climate changes 
induced by changes in atmospheric composition (increasing abundance of greenhouse gases and 
aerosols) and, consequently, cannot be ignored in evaluating the effects of cloud-radiation 
interactions on the climate. Identifying the sources of discrepancies will also help to identify the 
problems associated with the two approaches. Unfortunately, examining and isolating the sources 
of disagreement is complicated by the effects of different temporal and spatial sampling of the 
ISCCP and ERBE measurements. 

The existence of both narrowband and broadband channels on the ScaRaB radiometer 
provides a unique opportunity to compare the narrowband-based ISCCP and broadband-based 
ERBE approaches using coincident, collocated, and coangular radiance measurements and flux 
estimates. The narrowband channels of ScaRaB are roughly equivalent to the visible and infrared 
channels used in the ISCCP cloud retrievals, while the broadband channels are similar to those of 
the ERBE scanners. Thus, it is possible to use the ScaRaB data to retrieve the OLR using the 
ISCCP approach with the narrowband radiances and using the ERBE approach with the broadband 
radiances. Use of the ISCCP-like and ERBE-like OLRs from the ScaRaB measurements 
eliminates the uncertainties due to the temporal and spatial mismatches and allows us to focus on 
the differences between the ISCCP and ERBE approaches which is invaluable for the validation 
and improvement of both approaches and will significantly advance our understanding of ERB as 
a consequence. Particularly, simultaneous narrowband and broadband measurements allow us to 
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isolate the effects of angle model differences in the analyses'. Since our comparisons are restricted 
to the instantaneous pixel-level radiance and flux values, the differences in the space-time 
averaging algorithms used to obtain proper averages of the fluxes from temporally and spatially 
sparse samples are not considered in this study. There are other studies designed for this purpose 
(e.g., Xu 1997). Moreover, this particular study is only intended to identify the sources of 
disagreements in the ERBE and ISCCP approaches. The radiative effects of clouds using ISCCP 
and ERBE data are the subjects of other research efforts (e.g., Harrison et al. 1990; Hartmann and 
Doelling 1991; Ockert-Bell and Hartmann 1992; Hartmann et al. 1992; Chen et al. 2000a, b). 

A short description of the ScaRaB project and its ERBE-Iike inversion method is given in 
Section 2. The ISCCP approach to determine OLR using the ScaRaB narrowband radiances is 
described in Section 3. Section 4 relates the cloud cover categories from the broadband-based 
ERBE maximum likelihood estimation technique to the cloud types from the narrowband-based 
ISCCP cloud detection and retrieval algorithm. Section 5 examines the differences between the 
ISCCP-like and ERBE-like OLRs as a function of ISCCP cloud parameters and discusses the 
possible sources of disagreements. Section 6 summarizes the main findings and discusses the 
implications of this study. 


1 The unique advantages for addressing angular model related issues provided by ScaRaB have also been exploited in 
the study by Chang et al. (2000b), which focused on the TOA albedo and SW angular dependence models. 
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2. ScaRaB project and inversion method 

a. ScaRaB project 

The Scanner for Radiation Budget (ScaRaB) project is a joint project of France, Russia, 
and Germany intended to provide a continuation of the ERBE scanner mission (Kandel et al. 
1994). The first ScaRaB flight model (FM1) was launched on board the Russian Meteor-3/7 
weather satellite with circular polar orbit at an altitude of 1200 km and an inclination of 82.5° (not 
sun-synchronous). The orbit plane precessed relative to the sun-earth vector with a period of 
approximately 213 days, resulting in full coverage of the diurnal cycle in about 106 days. The 
ScaRaB instrument is a cross-track scanning radiometer with a swath of 97.82° (i.e., the maximum 
satellite nadir angle is 48.91°, yielding a maximum viewing zenith angle of 63.57° considering the 
curvature of the earth; Capderou 1998). The angular resolution of instantaneous field of view is 48 
x 48 mrad 2 , corresponding to a pixel size of 60 x 60 km 2 at nadir. The radiometer consists of two 
narrowband and two broadband channels with identical fields of view. The two narrowband 
channels are similar to the weather satellite imager channels used in the ISCCP analysis, one 
corresponding to visible portion of the solar spectrum (VIS: 0.55-0.65 pm), the other to the 
infrared atmospheric window (TRW: 10.5-12.5 pm). The broadband shortwave (SW: 0.2-4 pm) 
and total radiation (TW: 0.2-50 pm) channels of ScaRaB, from which the reflected shortwave and 
emitted longwave radiances are derived, are analogous to those of the ERBE scanner. Daytime 
radiation in the longwave band (LW: nominally 4-50 pm) is obtained by appropriated weighted 
subtraction of the SW signal from the TW signal. The radiometric accuracy of the radiances is 
estimated to be better than 1% in the LW and 2% in the SW domains. The data was received and 
processed from 24 February 1994 until 6 March 1995 with some interruptions; therefore, roughly 
1 1 months of data has been collected (Kandel et al. 1998). 

b. Inversion method 

As a follow-on of ERBE, the first version of the ScaRaB data processing was developed as 
an "ERBE-like version" to ensure consistency between ERBE and ScaRaB products. In other 
words, the ScaRaB algorithms are as close as possible to the ERBE algorithms except for some 
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parameters used in ERBE scene identification, which are missing in the published descriptions of 
ERBE (Yiollier et al. 1995) 2 , and the steps specific to the instrument characteristics, such as the 
spectral corrections. Like all scanning radiometers, ScaRaB measures radiances at specific angles. 
The filtered SW and LW radiances, after radiometric calibration, are first corrected for spectral 
filtering effects. Spectral correction is in general a function of scene type. For the ScaRaB LW 
radiance, which is the result of weighted subtraction of the SW from the TW channels, however, 
the filtering error is negligible and therefore no spectral correction is applied in the current 
ScaRaB inversion process (Viollier et al. 1995). The ERBE angular dependence models (ADMs), 
which account for the anisotropy of the radiation field, are then applied to the unfiltered SW and 
LW radiances to obtain the SW and LW fluxes at TOA (Suttles et al. 1988; Suttles et al. 1989). 

Since the angular dependence models depend on the type of surface and the amount of 
cloud present in the observed scene, a procedure is required to identify the scene type. The 
algorithm of scene identification is the maximum likelihood estimation (MLE; Wielicki and Green 
1989), which compares the unfiltered LW and SW radiances to a priori statistics developed from a 
classification of Nimbus-7 ERB data. The 12 scene classifications are the same as used by ERBE, 
based on a combination of five geotypes (ocean, land, snow, desert, coast) and four cloud 
categories: clear [0%-5% cloud cover], partly cloudy [5%-50% cloud cover], mostly cloudy [50%- 
95% cloud cover], overcast [95%-100% cloud cover]. Very roughly, the MLE expects a clear 
scene to be dark and hot, cloudy scenes to be bright and cold. Presumably, the larger pixel size of 
ScaRaB (60 km at nadir, versus 30-45 km for ERBE) should give rise to a different scene type 
distribution compared with ERBE. This effect, however, is mitigated because the a priori 
probabilities used in the ScaRaB scene identification were computed by analyzing the scene 
fractions of the ERBE S-9 archive (Viollier et al. 1995). Since the objective of the ScaRaB project 
is to determine the TOA ERB components on a scale appropriate to the study of cloud radiation 
interactions, an additional procedure is performed to obtain proper averages of the fluxes from 

2 Several improvements are being planned for the future version of ScaRaB data processing, including adapting 
ISCCP cloud detection algorithms to the coarse spatial resolution of the ScaRaB narrowband measurements for better 
scene/cloud identifications and thus improved clear-sky fluxes and cloud radiative effects (Briand et al. 1998) and 
applying different angular dependence models from recent studies for the conversion of radiances into fluxes (e.g, 
Stubenrauch et al. 1993). 
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temporally and spatially sparse samples (Brooks et al. 1986). However, as mentioned before, the 
problems associated with the spatial and temporal averaging algorithms are beyond the scope of 
this present study, since our results are restricted to the instantaneous pixel-level radiance and flux 
values. 

Similar to ERBE, three levels of ScaRaB data are produced by the data processing 
procedures. In particular, the second level (A2) is directly comparable to the ERBE S-8 product 
containing instantaneous unfiltered SW and LW radiances, estimated TOA fluxes, as well as scene 
identifications at the pixel level, and is the primary data source for this study (Kandel et al. 1998). 
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3. ISCCP-based approach for OLR determination 

The ScaRaB visible (VIS) and infrared window (IRW) channels are comparable to those 
used in the ISCCP analysis; thus it is a natural idea to apply the ISCCP algorithm to the ScaRaB 
VIS and IRW radiances. Similar to many cloud algorithms, the ISCCP algorithm consists of two 
basic phases: the cloud detection phase and the cloud properties derivation phase. Based on a 
series of statistical analyses of the radiance variations, the cloud detection algorithm first infers the 
clear-sky VIS and IRW radiances for each individual pixel. An individual pixel is then classified 
as cloudy, if either the VIS or IRW radiance differs from its corresponding clear-sky value by 
more than the detection thresholds, or clear if otherwise (Rossow and Garder 1993a). If the pixel is 
cloudy, then comparison of the observed radiances to those predicted by a radiative transfer model 
determines a cloud top temperature (T t ) from the IRW radiance and a visible optical thickness (T, 
daytime only) from the VIS radiance. If the pixel is clear, then surface temperature (T s ) and 
surface visible reflectance (R s , daytime only) are retrieved (Rossow et al. 1996). The effects of the 
atmosphere on the radiances are accounted for and cloud top pressure (P r ) is determined from T, 
using the atmospheric temperature and humidity profile from the Television Infrared Observation 
Satellite (TIROS) Operational Vertical Sounder (TOVS) product (Kidwell 1995). To be consistent 
with the cloud and surface retrieval, the same TOVS profile is also used in the radiative transfer 
model calculation in this study (unless otherwise indicated). 

From the ScaRaB VIS and IRW radiances, cloud and surface properties are retrieved in the 
same form as the ISCCP DX (pixel level). The ScaRaB-derived cloud and surface properties, 
called ScaRaB DX, are then combined with ancillary data and used as input to a radiative transfer 
model to calculate fluxes. The radiative transfer model is a refined version of the one described by 
Zhang et al. (1995) with similar input and output quantities. The improvements most relevant to 
this study are as follows. 1) The clouds are treated using an ice microphysical model if top 
temperature < 260 K, which is consistent with the ISCCP retrieval. 2) The spectral resolution of k- 
distribution method is increased from 25 to 33 k-intervals for thermal radiation (Lacis and Oinas 
1991). 3) The parameterization of water continuum absorption is improved using the theoretical 
calculations by Ma and Tipping (1991, 1992a, b). 4) Non-unit-emissivity is used for the surface 
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with correction of the retrieved surface skin temperature. 5) Better interpolation for surface air 
temperature is introduced. 

For the LW radiation, which is the focus of this study, the radiative transfer model uses a 
three-point quadrature (corresponding to emission angles of cos6- 1.0, 0.5, and 0.1) to integrate 
over the angular distribution of the radiation (Lacis and Oinas 1991). A by-product of this added 
complexity is that we have available from the radiation model an infrared window radiance, a 
spectrally integrated broadband radiance (both for nadir-viewing geometry and also for the more 
extreme emission angles of cosO — 0.5 and 0.1), and angle-integrated fluxes. Thus, the radiative 
transfer model can be used to simulate the narrowband and broadband ScaRaB measurements for 
direct comparison with the measured radiances, as well as allowing comparison of the resulting 
fluxes (Fig. 1). This aspect of the radiative transfer model makes the ISCCP approach well-suited 
to examine the effect that changing cloud properties have on the determination of ERB fluxes. 

Although all 1 1 months of ScaRaB FM1 narrowband radiance data have been processed 
through the ISCCP cloud algorithm, our analyses are performed on only two months of ScaRaB 
DX data since, typically, one month of ScaRaB DX data consists of several million observations 
and hence are enough to produce statistically significant results. However, since the satellite 
retrieval and flux calculation/conversion also depends on the relative geometry among sun, Earth 
target, and satellite, we choose two specific months, namely May and July 1994, to provide a 
representative sample of sun-target-satellite geometries. In May 1994, the daylight Equator 
Crossing Time (ETC) of the satellite was around early afternoon, drifting from about 1500 local 
time (LT) at the beginning of the month to about 1200 LT at the end of the month. Because of the 
precession of the orbit, the daylight ETC changed to early morning in July 1994, varying from 
0900 to 0700 LT approximately. Only the daytime observations are used in this study, since both 
approaches are more accurate for the daytime than for the nighttime observations due to the fact 
that, lacking visible information, both the ISCCP cloud detection and the ERBE MLE scene 
identification algorithms degrade into single threshold tests for the nighttime observations. 
Besides, in the ISCCP retrieval, cloud optical thickness is determined from visible radiance 
measurements, which are only available during daytime, and the cloud top temperature is more 
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accurately determined as a function of optical thickness (Rossow et al. 1996; Rossow and Schiffer 
1999). 

One of the problems encountered when applying the TSCCP algorithm to the ScaRaB 
measurements arises from the coarse spatial resolution of ScaRaB measurements (60 km at nadir), 
since the ISCCP algorithm is optimized for use with pixel size of 4-7 km. Differences in the pixel 
size primarily affect the retrieval in two ways. 1) The spatial and temporal variability of the 
surface characteristics decreases as the pixel size increases. This will cause some cloud- 
contaminated ScaRaB pixels to be declared clear by the original ISCCP cloud detection algorithm, 
where the thresholds in the space/time variability test are selected for use with higher spatial and 
temporal variations of radiances associated with the finer spatial resolutions (Briand et al. 1998). 
2) ISCCP assumes that the cloud optical properties are uniform at pixel scales; hence, the cloud 
cover of pixels is assumed to be either zero or one. This assumption becomes less acceptable at the 
lower ScaRaB spatial resolution. To separate the clear (cloudy) ScaRaB DX pixels with minimal 
cloud (clear-sky) contamination from the clear (cloudy) pixels with possible cloud (clear-sky) 
contamination, the ScaRaB DX pixels are compared with the actual ISCCP DX pixels matched in 
both space and time as closely as possible. The comparison divides the matched ScaRaB DX 
pixels into two subsets, called "unmixed" and "mixed". The clear (cloudy) ScaRaB DX pixel is 
classified as "unmixed" only if its four nearest (within a radius of 60 km) 3 neighboring ISCCP DX 
pixels are all labeled as clear (cloudy) and is classified as "mixed" otherwise. For the "unmixed" 
pixels, we assume that when four nearby ISCCP DX pixels are clear (cloudy), the whole region of 
ScaRaB DX pixel is entirely clear (cloudy). Conversely, the "mixed" pixels are more likely to be 
partially cloud covered. On average, about 65% of the matched clear (determined by ISCCP 
analysis) ScaRaB DX pixels are labeled as "mixed", indicating that a quite large portion of the 
ISCCP-declared clear ScaRaB pixels may actually contain small amount of clouds to compensate 
for the effects of the lower spatial resolution, as pointed out by Briand et al. (1998). A large 

3 Since the ScaRaB pixel size grows to over 200 km at 60° view, limb-viewed ScaRaB pixels should be compared 
with significantly more neighboring ISCCP DX pixels in theory. However, if more ISCCP DX pixels were used, the 
number of “unmixed” pixels would have decreased significantly. To retain enough “unmixed” pixels for analysis, four 
nearest ISCCP DX pixels are used regardless of the ScaRaB viewing angle. As a result, a portion of limb-viewing 
“unmixed” ScaRaB pixels are in reality “mixed”. The consequences are considered later in Section 5. 
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majority (about 75%) of the ISCCP-declared cloudy ScaRaB pixels are labeled as "unmixed" 
rather than "mixed", especially for the cloud types that are optically thick or of large size such as 
deep convective or cirrostratus clouds (see Fig. 2 for the ISCCP radiometric cloud classification). 
The only exception arises for the ScaRaB pixels determined as cumulus by the ISCCP algorithm, 
where more than half (roughly 58%) of them are labeled as "mixed", consistent with the facts that 
1) cumulus clouds tend to be small-sized and broken; and 2) partially-filled cloudy pixels are more 
likely classified as cumulus rather than any other cloud type in the ISCCP analysis. Although 
completely clear or overcast does not imply that the optical properties of surface or cloud are 
uniform over areas of tens of kilometers, the assumption of optical homogeneity is certainly a 
better one for the "unmixed" pixels than for their "mixed" counterparts. It is therefore interesting 
to contrast the results from the "unmixed" pixels with those from the "mixed" ones and examine to 
what extent the results are affected when this assumption is severely violated, as will be done in 
section 5. 
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4. Relationship of MLE cloud cover categories to ISCCP cloud types 

In the ERBE MLE, the full range of observed cloud height and optical thickness variations 
are only implicitly included in the derived cloud cover categories (Wielicki and Green 1989). 
ISCCP cloud detection and classification, on the other hand, attributes all variations in the 
narrowband radiances to changes in cloud top temperature and visible optical thickness, assuming 
clouds cover individual pixels completely. Examination of the MLE-derived cloud cover 
categories within the context of ISCCP-derived cloud types, defined by the cloud top pressure and 
visible optical thickness (Fig. 2), will not only allow us to better understand the MLE-derived 
cloud cover, but also help us to explain the differences between the ERBE- and the ISCCP- 
approach radiative fluxes. 

Tables 1 and 2, respectively, give the May 1994 monthly global results of the distribution 
of the ISCCP-derived scene types in each of the four MLE-derived cloud cover categories and the 
distribution of the four MLE-derived cloud cover categories in each of the ISCCP-derived scene 
types. The ISCCP-derived global mean ScaRaB cloud amount (62.9%), defined as the percentage 
of cloudy pixels, is about 5% lower than the 8-yr (1986-93) global mean cloud amount of the 
actual ISCCP D-series data (67.6%, Rossow and Schiffer 1999). Since the ScaRaB spatial 
resolution is much lower than the original ISCCP resolution, this result strongly indicates that the 
ISCCP algorithm is quite stable with spatial resolution variations, as also observed by Briand et al. 
(1998). Compared with the ISCCP-derived ScaRaB clear-sky amount (37.1%), the MLE-derived 
ScaRaB clear-sky amount is about 11% less (Table 2). This is not unexpected since, on the one 
hand, the original ISCCP algorithm, when applied to the coarser ScaRaB measurements, tends to 
overestimate the clear-sky amount by misidentifying some cloud-contaminated pixels as clear 
(Briand et al. 1998) and on the other hand, the MLE clear-sky amount is biased low because the 
definition of clear sky in the ERBE analysis is overly restrictive (Collins and Inamdar 1995). For 
the ScaRaB pixels declared as clear by the ISCCP cloud detection algorithm, more than half of 
them are also labeled as clear by the MLE scene identification scheme, and less than 10% of them 
are labeled as overcast by MLE. 
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However, the most distinctive feature shown in Table I is the progressively higher MLE- 
derived cloud cover for the ISCCP-derived cloudy pixels with larger optical thickness or lower top 
pressure values because of the stronger solar or thermal contrast between these clouds and clear 
sky, consistent with Diekmann and Smith (1989). This dependence implies that the MLE-derived 
cloud cover is likely to be biased low for low thin clouds and biased high for high opaque clouds. 
In other words, the MLE-derived cloud cover will best approximate true cloud cover for clouds 
with average height and optical thickness, since the a priori data used by the MLE is an average 
over all observed cloud heights and optical thicknesses, as pointed out by Wielicki and Green 
(1989). Specifically, for the ScaRaB pixels classified as cumulus by ISCCP, characterized by their 
relatively weak radiative contrast from clear sky, they are most likely classified as partly cloudy 
by MLE. For the ScaRaB pixels classified as stratocumulus, stratus, altocumulus, or altostratus by 
ISCCP, the MLE scene identification most likely declares them as mostly cloudy. This result is 
consistent with Chang et al. (2000a), who also find that nearly all of the uniform low-level, single- 
layered, marine overcast clouds (based on the spatial coherence analysis applied to AVHRR 
measurements) were classified as mostly cloudy by ERBE MLE. Finally, for the ScaRaB pixels 
classified as cirrus, nimbostratus, cirrostratus, or deep convective cloud by ISCCP, most likely 
they are declared to be overcast in the MLE, showing that the MLE method expects the overcast 
pixels to be both cold and bright by design. Overall, for the ScaRaB pixels declared as cloudy by 
ISCCP, the MLE method most likely classifies them into mostly cloudy scenes. 

In the MLE classification, for May 1994, the ScaRaB pixels are roughly evenly distributed 
among the four MLE cloud cover categories, where the difference between highest and lowest 
frequency of occurrence (FOC) is less than 9% (Table 2). Among the ScaRaB pixels classified as 
clear by MLE, more than 80% of them are also declared clear by ISCCP. For the ScaRaB pixels 
classified as partly cloudy by MLE, more than 90% of them are declared either clear sky, low 
clouds with thin or medium optical thickness (cumulus and stratocumulus), or thin middle-level 
clouds (altostratus) by ISCCP. For the MLE-derived mostly cloudy and overcast ScaRaB pixels, 
the ISCCP algorithm tends to classify them as clouds with even larger optical thickness and/or 
lower top temperature. Again, this table shows that the ERBE MLE-derived cloud cover implicitly 
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includes information on cloud height and optical thickness, producing mixtures of r-variation and 
P, -variation. Since the ERBE ADMs are empirically constructed by sorting the Nimbus-7 ERB 
measurements into latitude bins, angular bins and cloud cover categories, it is anticipated that for 
each particular cloud cover category, the associated angular model tends to be biased to the 
angular behavior of the ISCCP scene types with the highest FOCs. The most likely scene types, 
however, may differ from the ones in this study, since the cloud detection algorithm used to 
construct the ERBE ADMs was the new cloud-ERB (NCLE) algorithm (Stowe et al. 1988) and 
significant differences exist in cloud fraction between the MLE and NCLE results (Suttles et al. 
1992). It is worth noting that considerable differences in cloud fraction exist between the ISCCP 
and NCLE results, too (Rossow et al. 1993). 

A closer look at the two classifications reveals that the relation between ISCCP- and MLE- 
derived scene types varies with the satellite viewing zenith angle (VZA). Figure 3, showing the 
variation of MLE-derived scene type frequency of occurrence (FOC) as a function of VZA 
(grouped by ISCCP-determined scene types) for May 1994, reveals that the MLE scene type FOC 
varies considerably with the satellite viewing angle for some particular ISCCP-derived cloud 
types. Most notably, for the ScaRaB pixels identified as cirrus by ISCCP, the FOC of overcast in 
the MLE increases from about 20% for VZA = 0°-15° to over 40% for VZA > 51°. This 
considerable variation may result from the viewing angle dependence of MLE scene identification 
scheme (Ye and Coakley 1996a, b) or from that of ISCCP cloud detection and retrieval algorithm 
(Rossow and Garder 1993b), but more likely from both. The viewing angle dependence of MLE 
implies that the two steps used in the ERBE approach to derive the radiative fluxes, that is, first 
identifying within each scanner field of view the scene type, then applying an anisotropic factor 
associated with that scene type to convert radiances into fluxes, are intricately linked to each other 
and should not be viewed separately. When the entire population of pixels is included, however, 
no viewing angle dependence is apparent in the MLE scene type FOC (except for a very slight 
decrease in clear scene FOC as VZA increases) regardless of the month considered. In reality, 
since the ScaRaB field of view size grows with the viewing zenith angle, from ~60 km at nadir to 
over 200 km at limb, one would expect that FOC increases for partly and mostly cloudy scenes 
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and decreases for clear and overcast scenes. This is not the case in Figure 3, which suggests that 
the scene type names used by the MLE classification are rather inexact and misleading, and hence 
should not be interpreted literally. Similar results are also obtained for the near-terminator month, 
July 1994 (not shown). 
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5. Comparison between the two approaches 

As discussed earlier, both the ISCCP retrieval algorithm and radiative transfer model for the 
flux calculation assume optical homogeneity in the underlying pixels. Undoubtedly, this 
assumption becomes increasingly a problem as the pixel size increases. To alleviate this problem 
associated with the relatively coarse spatial resolution of ScaRaB measurements, we will first 
restrict the results to the “unmixed” scenes. Then the results from the “mixed” scenes will be 
briefly discussed to determine how the results are affected when the optical homogeneity 
assumptions are severely violated. Apart from the above mentioned optical homogeneity issue, our 
comparisons are complicated further by the following two differences between the ERBE and 
ISCCP approaches. 1) The TOA is chosen as 30 km in the ERBE analysis, whereas the model 
calculations set TOA as 100 km. Using 30 km instead of 100 km as the definition of TOA in the 
radiative transfer model produces a fraction of a percent difference in the calculated radiances or 
fluxes. 2) The synthetic ScaRaB LW channel (the result of weighted subtraction of the SW from 
the TW channel) has a nominal spectral range from 4 to 50 pm, whereas the model calculated LW 
radiance or flux represents the range from 5 to 200 pm. This difference may cause a low bias of 
as much as 1% in the ScaRaB LW measurements. 

a. “Unmixed" pixels 

Figs. 4a and 4b show the mean differences (over ocean, “unmixed” pixels only) between 
the measured and model-simulated TOA LW radiances, between the converted TOA LW fluxes 
using measured LW radiances (results from two different ADMs, called ERBE and SDK ADMs 
— see later text for details) and model calculated fluxes for May and July 1994, respectively. 
Similar results are obtained over land (not shown). Results from clear and cloudy, from nadir- and 
limb-viewing pixels (VZA = 0° and 60°) are separately presented. Cloudy pixels are further 
grouped into 15 classes according to the ISCCP retrieved cloud top pressure (separated at 680 mb 
and 440 mb) and optical thickness (divided at 1, 3.6, 9.4, and 23). Except for the clear and low 
cloud cases, where the differences exhibit a simple latitudinal dependence (maximum in the 
Tropics), the differences are generally independent of the geographical location, and hence only 
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the global mean results are discussed (likewise the numbers quoted in the text, unless otherwise 
noted). Figs. 5a and 5b are similar to Figs. 4a and 4b, except that only the pixels where the 
measured and model-simulated LW radiances agree to within 0.5 Wm' 2 sr 1 are included, so that 
the LW flux difference can originate only from the difference between the ADMs used by the 
ERBE-like approaches and the radiative transfer model’s treatment of the angle integration. 

i. Clear sky and low clouds 

Even for the “unmixed” clear ScaRaB DX pixels, a significant portion of them (26.2% for 
May 1994, and 44.0% for July 1994) are identified as somewhat cloudy by the ERBE MLE 
method. To minimize the uncertainties in clear-sky results caused by the possible cloud 
contamination, only the clear “unmixed” pixels where both ISCCP and MLE agree are included in 
the comparison. Averaged over the globe, the clear-sky TOA LW fluxes estimated from the 
observed LW radiances and ERBE ADMs always exceed the model calculations, ranging from ~3 
Wm" 2 for limb-viewing (VZA = 60°) land pixels of July 1994 (not shown), to ~14 Wm 2 for nadir- 
viewing ocean pixels of May 1994. Note that the high bias in the ERBE analysis shown here 
cannot be explained by the conventional wisdom that the ERBE clear sky TOA LW fluxes are 
biased high, especially in the Tropics, due to the overly restrictive clear sky definition in the 
ERBE analysis (cf. Hartmann and Doelling 1991; Collins et al. 1995) for the following two 
reasons: 1) the current flux comparison is made on the same “unmixed” clear pixels where ISCCP 
and MLE agree; and 2) these “unmixed” clear pixels, by definition, tend to be found in the midst 
of vast regions that are clear, and thus are highly unlikely to be cloud-contaminated. In light of the 
fact that ERBE ADMs do not explicitly take into account the attenuation of LW radiation by the 
atmosphere above the surface or cloud top, Stubenrauch et al. (1993) parameterized the longwave 
ADMs using combined narrowband IRW and broadband LW observations to specifically include 
the effects of cloud height variations (referred to as SDK ADMs in this paper). Replacing the 
original ERBE ADMs with the SDK ADMs only slightly reduces the discrepancy in the clear sky 
TOA LW fluxes. The comparison between the observed and model-simulated LW radiances 
shows that the flux differences largely arise from the disagreements in the broadband radiances. 
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This is further confirmed in Figs. 5a and 5b. For the clear pixels where the measured and 
simulated radiances agree, the flux differences between radiance converted (using either ERBE or 
SDK ADMs) and model calculated fluxes are very small (in the worst case 2.5 WnT 2 ), indicating 
that these two ADMs (especially the ERBE ADMs) and the radiative transfer model’s 
representation of longwave emission anisotropy are fairly similar for clear conditions. 

A detailed investigation of this clear sky difference shows that there is a positive 
correlation between the LW radiance differences (observation minus calculation) and the column 
precipitable water amount (Fig. 6). Considering the sensitivity of the radiative transfer model to 
the input atmospheric temperature and humidity profiles, we repeated the calculations using the 
atmospheric profile from the TOVS Pathfinder path-B dataset (hereafter 31 profile, from the name 
of the retrieval algorithm used; Scott et al. 1999). The use of 31 profile brings the simulated LW 
radiances closer, in a global mean sense, to the observations by 1-2 Wm 2 sr ', leading to a 3-6 Wm 
2 reduction in the flux difference; however, a stronger correlation between LW radiance 
differences and the column precipitable water amount is observed, regardless of the month and 
viewing zenith angle, implying a larger difference between the ADM-derived and model 
calculated latitudinal gradients of radiative fluxes using the 31 than the TOVS profile. Note that 
this result does not necessarily suggest that the water vapor itself (either the errors in the input 
humidity data or the treatment of water vapor in the model) is the source of the discrepancy, 
because other properties of the atmosphere that affect longwave radiation (such as temperature) 
are highly correlated with water vapor. 

We have repeated the radiative transfer calculation by adjusting several model inputs to 
which the TOA longwave radiation is most sensitive, such as increasing the atmospheric or 
surface temperature, or reducing the atmospheric humidity within their assessed uncertainty ranges 
(cf. Zhang et al. 1995). Though each of the adjustments brings the model calculation closer to the 
observation, none of them is, by itself, sufficient to explain the large discrepancy between model 
calculation and observation for the nadir-viewing ocean pixels of May 1994. It is likely that the 
discrepancy arises as a combined result of uncertainties in the model inputs and/or errors in the 
radiative transfer model. However, since very different amount of changes to the model inputs are 
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needed to get agreement between model calculations and observations in May and July, the 
discrepancy in LW radiances/fluxes may be due more to uncertainties in the model inputs than to 
errors in the radiative transfer model. Those uncertainties/errors need to be investigated and 
quantified with further research. 

Similar to the clear-sky results, but at a reduced magnitude, both the observed LW 
radiances and the radiance-converted fluxes exceed the model calculations for the ScaRaB pixels 
identified as containing low clouds (P, > 680 mb) by ISCCP (regardless of the MLE-derived cloud 
cover), due mainly to the clear-sky effect. Figure 5 shows that the TOA LW flux differences for 
low clouds stemming from the different angle treatments are generally very small as well, despite 
a weak but systematic dependence on optical thickness and viewing zenith angle. Similar but 
stronger dependencies are also observed for middle- and high-level clouds, which will be 
discussed next. 

ii. High- and middle-level clouds 

As shown in Figs. 4a and 4b, for the ScaRaB pixels classified as high clouds (P r < 440 mb) 
by ISCCP, the model-simulated LW radiances can be either less or greater than the observed ones, 
depending on the retrieved optical thickness and viewing zenith angle. When the clouds are very 
optically thin, the model in general underestimates the TOA LW radiance. At nadir-viewing 
geometry, the difference is about 4 Wm' 2 sr' 1 . Analysis shows that the bias primarily arises from 
the visible adjustment procedure in the ISCCP retrieval. ISCCP first retrieves the cloud top 
temperature assuming that the cloud is opaque to infrared radiation, then adjusts the retrieved 
cloud top temperature if the infrared transmission of the cloud estimated from the retrieved visible 
optical thickness is greater than 0.5% 4 . This procedure encounters difficulties with very thin 
clouds because the errors in the optical thickness, surface temperature, and the ratio used to relate 
the optical thickness value in the VIS to that in the IRW will be greatly amplified when the 


4 A small coding error in the ISCCP visible adjustment step that causes a slight overestimation of cloud top 
temperature and pressure for clouds with visible optical thickness between 2 and 6 has been found. This error has been 
corrected for the ScaRaB DX data used in this study and for the regular ISCCP D data for 1994 and later. However, 
ISCCP D data for July 1983 through 1993 still contain this error (see appendix). 
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infrared transmission is too large. Thus, for the clouds with xlcosQ < 1.2, ISCCP tends to place 
them at the tropopause and reset the optical thickness to a proper value so that the infrared window 
(but not the VIS) radiance is preserved (Rossow et al. 1996). This treatment inevitably leads to a 
cold bias in the calculated broadband longwave radiance, due to the different cloud-radiation 
interactions inside and outside the infrared window wavelength range: even though the 
combination of T ( and r reported by ISCCP matches the radiation in the window, if the clouds are 
too high they reduce the radiation outside the window too much. 

As for the high clouds with medium or large optical thickness, the model-simulated LW 
radiances agree rather well with the observations for nadir-viewing pixels, but always exceed the 
observations for pixels with VZA = 60°. Noticing that the agreement is worse for the “mixed ’ 
pixels (cf. Fig. 8), we speculate that the larger discrepancy at 60° view is caused, at least in part, 
by forcing the ISCCP analysis (optimized for use with pixel size of 4-7 km) to operate on the 
excessively large limb-viewing pixels that tend to have a larger degree of sub-pixel cloud 
heterogeneity and brokeness (even for the “unmixed” pixels). In fact, we have repeated the 
analysis by comparing the ScaRaB DX pixels with the 16 nearest ISCCP DX pixels (within a 
radius of 200 km, roughly the size of limb- viewing ScaRaB pixel) and found that more than half 
of the original “unmixed” ScaRaB DX pixels are no longer “unmixed”. For those ScaRaB pixels 
that remain “unmixed”, the standard deviation of cloud top temperatures over the 16 neighboring 
ISCCP DX pixels ( crj ,) is computed and used as the indicator of sub-pixel cloud heterogeneity 
within each individual ScaRaB DX pixel (cf. Stubenrauch et al. 1993). We found that both the 
mean and standard deviation of the discrepancy between the model-simulated and observed LW 
radiances at 60° view decrease as err, gets smaller (Fig. 7). 

Overall, the model-simulated radiances agree better with the measurements for the middle- 
level clouds (440 mb < P, < 680 mb). However, this agreement is more likely due to the 
cancellation of two opposing effects, namely the “clear sky” bias (mostly affecting clear sky and 
low clouds) and “sub-pixel cloud heterogeneity/brokeness” effects (mostly affecting limb-viewed 
high clouds), with both effects operating at a smaller magnitude. As a result, the vertical 
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distribution of the longwave cooling and surface radiation might not be correctly simulated, even 
if the TOA radiation is correctly calculated for the midlevel clouds. 

i ii. ERBE and SDK ADMs and the model’s treatment of angle integration 

The fact that the differences in ADM-converted and model calculated flux only 
approximately follow the radiance differences suggests that part of the flux differences are caused 
by the difference between the ADMs used and the radiative transfer model’s treatment of angle 
integration. Although the relative contribution of the different angle treatments to the flux 
differences is difficult to determine exactly, Figs. 5a and 5b do provide some insight. 

We first consider the ERBE ADMs. It is observed in both Figs. 5a and 5b that for the nadir 
view, when the observed and the model-simulated radiances agree, the ERBE ADM converted 
fluxes are higher (lower) than the model calculated ones for clouds with small (large) optical 
thickness, regardless of the cloud height. The opposite is true for the 60° view. This indicates that, 
compared with the radiative transfer model, the limb-darkening in the ERBE angular models is too 
weak for the optically thin clouds (r< 3.6), but too strong for the optically thick clouds (r> 9.4). 
When all clouds are included, however, the limb-darkening in the ERBE angular models closely 
resemble that implied by the radiative transfer model (not shown). In reality, for clouds at the 
same height, optically thin clouds usually are more limb-darkened than the optically thick clouds. 
For instance, due to the seed dependence of radiation path length through the cloud layer, optically 
thin cirrus can be nearly transparent to LW radiation at near nadir viewing angles, while the same 
cloud can be nearly opaque to LW radiation at near-limb viewing angles, leading to the limb- 
darkening even stronger than that for clear-sky conditions (Wielicki and Green 1989). On the 
other hand, the optically thick clouds are usually more Lambertian compared with clear scenes. 
These results suggest again that the ERBE ADMs represent the anisotropy of the LW radiation 
averaged over all cloud types but not for particular cloud types. This is not surprising considering 
that the ERBE ADMs are based on a large set of Nimbus-7 ERB observations without explicit 
distinction of different cloud optical properties. As pointed out by Stubenrauch et al. (1993), if the 
dependence of angular models on the cloud properties was not taken into account in their 


21 


construction, it must be averaged out in the final products. In fact, this is exactly what happened, 
since the longwave ERBE ADMs are all quite similar with angle variations of only 1% to 2% on 
average among the four cloud conditions (Suttles et al. 1989). Also notice from Figs. 5a and 5b 
that the difference between the ERBE ADMs and the angle treatment of the radiative transfer 
model increases with cloud height. As a result, the differences in the estimated fluxes caused by 
the different angle treatments between the ERBE and ISCCP approaches are not the same for the 
low as for high clouds. This is very important since the radiative flux gradient, one of the ultimate 
forcings driving the atmospheric circulation, between the climate regimes characterized by 
persistent but different cloud types is not equally represented by these two approaches. Yet, this 
difference in the gradient cannot be completely eliminated in the final space- and time-averaged 
products, as will be discussed in the next section. 

Figs. 5a and 5b show that for the optically thick clouds (r> 9.4), the SDK ADMs and the 
radiative transfer model’s angle treatment are very similar. This is anticipated since the SDK 
ADMs are developed on the basis of radiative transfer calculations assuming a horizontally 
homogeneous surface or cloud layer of unit emissivity (Stubenrauch et al. 1993). Interestingly, for 
the optically thin clouds, the SDK ADMs appear to coincidentally behave more like the ERBE 
ADMs. 

b. ‘‘Mixed" pixels 

Figs. 8 and 9 are the same as Figs. 4a and 5a, respectively, except they show results for the 
“mixed” pixels. For high-level broken cloud fields, the agreement between model-simulated and 
observed LW radiances deteriorates. However, this effect is rather small on the global scale 
considering the relatively infrequent occurrence of high-level broken cloudiness: high clouds, 
usually large in size, tend to be continuous over a large area. In fact, the majority of the broken 
cloudiness is confined to the low level, small cumulus clouds in the trade wind regimes. For mid- 
and low-level clouds, the patterns shown in Fig. 8 are very similar to those in Fig. 4a, indicating 
that the uncertainties in the estimated TOA longwave radiation budget caused by broken 
cloudiness are quite small, even at the coarse spatial resolution of ScaRaB. This can be explained 
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by the cancellation of two offsetting effects. In the presence of broken cloudiness, the radiance 
used to retrieve the cloud is an averaged value from the actual clouds and surface, hence, the 
retrieved cloud assumed to cover the entire pixel is on average lower and optically thinner than the 
actual clouds. The larger cloud cover tends to decrease the longwave radiation whereas the higher 
cloud top temperature and smaller optical thickness tend to increase the longwave radiation. In the 
low cloud cases, the effect of broken cloudiness on the TOA longwave radiation budget is even 
smaller because the TOA longwave radiation is largely determined by the atmosphere above the 
clouds. 
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6. Summary and discussions 

Two conceptually different approaches (broadband-based ERBE and narrowband-based 
ISCCP approaches), used to derive the TOA longwave radiative fluxes, are compared using the 
ScaRaB simultaneous narrowband and broadband measurements (cf. Fig. 1). To evaluate the 
ERBE scene identification algorithm and facilitate the flux comparisons, the ERBE MLE-derived 
scene types are first compared with the ISCCP-derived cloud types on a pixel-by-pixel basis. The 
comparison shows a progressively higher MLE-derived cloud cover for the ISCCP-derived cloudy 
pixels with larger optical thickness or colder top temperature, indicating the cloud properties are 
implicitly included in the derived cloud cover categories (cf. Wielicki and Green 1989). Moreover, 
the comparison also confirms a view-zenith angle dependence of the MLE scene identification 
scheme (cf. Ye and Coakley 1996a, b), which could mean that the ERBE ADM converted radiative 
fluxes can be erroneous even if the ADMs themselves are correct, because they are not selected 
correctly at all viewing angles. But in reality, the errors in the scene identification scheme are 
more likely offset to some extent by the errors in the ADMs if the same scene identification 
scheme is used in the development and application of ADMs (Chang et al. 2000a). Thus, the scene 
identification and radiance-to-flux conversion steps in the ERBE processing are intricately linked 
together and should not be viewed separately. 

For the clear sky, the model underestimates the observed TOA LW radiances on the global 
average. Though the bias is found to be correlated on the column precipitable water amount, the 
exact source of the bias remains undetermined because of the possible dependencies on other 
variables, most notably temperature, that are correlated with water vapor variations. A similar but 
smaller bias is also observed for the low clouds, primarily due to the clear sky effect since the 
TOA longwave radiation is largely determined by the atmosphere above when the clouds are low- 
level. Further studies are needed to address this problem. For high clouds, depending on the 
optical thickness, the model-simulated LW radiances can be either less or greater than the 
observed ones. When the clouds are very optically thin, the model underestimates the TOA LW 
radiance as a result of the fact that these clouds tend to be placed too high during the visible 
adjustment step of ISCCP retrieval. In general, the modeled TOA LW radiances agree better with 
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observations for the nadir-viewing pixels than for the pixels viewed from limb, since the 
excessively large limb-viewing pixels are more likely to be heterogeneous or partially cloud-filled. 

Compared with the radiative transfer model, the limb-darkening in the ERBE angular 
models appears to be too weak for the optically thin clouds, but too strong for the optically thick 
clouds, and the difference between the ERBE ADMs and the angle treatment of the radiative 
transfer model increases with cloud height. Limitations in the radiative transfer model 
calculations, however, prevent definitive judgement on the ERBE angular models for the 
following reasons: 1) The radiative transfer model’s angle treatment is sensitive to the 
uncertainties in the cloud vertical structure (e.g., multi-layer clouds); and 2) the ERBE ADMs are 
generated using the complete range of observed cloud conditions, not just clear and overcast, 
whereas the plane-parallel radiative transfer model used in this study inevitably assumes optical 
homogeneity on the underlying pixels, which tends to underestimate the limb-darkening for the 
broken three-dimensional cloud fields. Nevertheless, this study clearly demonstrates that obtaining 
more accurate instantaneous longwave flux estimations from ERBE approach would require 
additional cloud classes based on cloud height and optical thickness as is planned for the advanced 
CERES analysis (Wielicki et al. 1996). 

One of the rationales behind the ERBE approach is that the viewing-zenith angle 
dependence of derived longwave flux caused by errors in the angular models will be significantly 
reduced in the final space- and time-averaged fluxes. This error reduction results from the fact that 
the ERBE angular models are normalized by the hemispherical flux based on the integral 
relationship between radiance and flux (Suttles et al., 1989). Over a long enough time periods, a 
satellite will view a given region of the earth from a wide range of viewing zenith angles, in effect 
averaging over the viewing zenith angle dependence of the derived longwave fluxes. However, the 
errors cannot be completely eliminated 1) because of the inability of a satellite platform to provide 
a uniform sampling in angle (Suttles et al. 1992), 2) because a particular radiometer may only scan 
over a limited range of viewing zenith angles (for example, the maximum viewing zenith angle of 
ScaRaB scanning radiometer is 63.57°) or the measurements are not used beyond a cutoff nadir 
angle since the measurements become increasingly difficult to interpret in terms of the radiant 
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exitance (Smith et al. 1986), and 3) because the possible convolution between the scene 
identification and radiance-to-flux conversion steps in the ERBE processing suggests that 
normalizing the angular model of an individual scene type independently may not be sufficient. 
Together with the fact that the magnitudes of the errors are cloud type dependent, these residual 
errors will certainly introduce systematic regional errors in the estimated mean and temporal 
variations of TOA LW radiation field. For example, as one of the better known interannual 
variations in climate, ENSO is characterized by large changes in cloud types in consequence of 
changes in the general circulation pattern of the tropical atmosphere. Unless these systematic, 
cloud type dependent errors are eliminated, one cannot accurately quantify the spatial and 
temporal variations in the TOA LW radiative fluxes associated with an ENSO event. Likewise, the 
ERBE-estimated TOA LW radiative flux gradient across the equatorial Pacific ocean needs to be 
revisited, considering the persistent and pronounced differences in cloud types over the ascending 
and descending branches of Walker circulation. 
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APPENDIX 


Error in the visible adjustment step of the ISCCP cloud retrieval 

For each cloudy pixel, ISCCP first retrieves the cloud top temperature assuming that the 
cloud is opaque to infrared radiation, then adjusts the retrieved cloud top temperature and visible 
optical thickness through iteration if the infrared transmission of the cloud is greater than 0.5% 
(Rossow et al. 1996). Unfortunately, a small coding error in the ISCCP cloud retrieval program 
that occasionally causes this iteration to stop prematurely has been found (for the May 1994 
ScaRaB DX, this coding error affects about 13% of daytime cloudy pixels). Even though the 
retrieval program was modified to rectify this error to produce the ScaRaB DX used in this study 
(as well as for the regular ISCCP data product beyond 1993), we document the error and assess to 
what extent it affects the retrieved cloud properties and calculated TOA radiative fluxes. 

Fig. Al shows the zonal mean difference between the ScaRaB DX cloud top pressure 
retrieved using the revised program and the program with the coding error in the visible 
adjustment step. By presenting the results from different optical thickness intervals separately, it is 
clear that this error is most noticeable for the cloudy pixels with medium optical thickness values 
(approximately from 2 to 6) and usually gives rise to a cloud top temperature/pressure value that is 
a bit too large (about 60 mb in P, in the worst case). Similar results were also obtained for one 
month (July 1994) NOAA DX data (not shown). It is also observed that this error decreases as 
cloud top temperature/pressure increases (not shown). Correction of this error has the effect of 
moving the cloudy pixels a little bit higher in altitude while also making them slightly optically 
thinner, resulting in only a very small shift in the relative abundance of the standard ISCCP cloud 
types. Analysis shows that the small shift occurs mostly from altostratus to cirrus clouds with the 
changes in either cloud type amounting to less than 1% (absolute value). 

Fig. A2 examines the resulting difference in the model calculated TOA LW radiative 
fluxes using one month (May 1994) ScaRaB DX data (revised minus flawed). This figure shows 
that the error leads to an overestimation of several Wm' 2 in the zonal mean TOA LW fluxes for 
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clouds with medium optical thickness values, consistent with the bias in the retrieved cloud top 
pressure. When all pixels (cloudy and clear) are included, the global mean TOA LW flux bias 
owing to this error is only 0.6 Wm' 2 , indicating that the overall impact of this error on the TOA 
radiation budget calculation is quite small. 
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Table captions: 


Table 1. Distribution of the ISCCP-derived scene types in each of the MLE-derived cloud cover 
categories, as well as the frequency of occurrence of each ISCCP-derived scene type for all the 
ScaRaB pixels spanning May 1994. The highest values are shown in bold. 

Table 2. Distribution of the MLE-derived cloud cover categories in each of the ISCCP-derived 
scene types, as well as the frequency of occurrence of each MLE-derived cloud cover category for 
all the ScaRaB pixels spanning May 1994. The highest values are shown in bold. 
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Figure captions: 


Figure 1. Schematic diagram of the basic approach used in the present study. 

Figure 2. Definitions of the nine different cloud types by ISCCP in terms of cloud-top pressure 
and visible optical thickness. The names used for convenience are meant to suggest qualitative 
relationships with the classic morphological cloud types, but they should not be interpreted to be 
quantitatively correct in every instance (cf. Lau and Crane 1995). 

Figure 3. Variation of MLE-derived cloud cover category (c: clear; p: partly cloudy; m: mostly 
cloudy; o: overcast) frequency of occurrence (FOC) as a function of viewing zenith angle for each 
of the nine ISCCP-derived cloud types, clear pixels, all cloudy pixels, and all pixels (data month: 
May 1994) 

Figure 4a. Mean differences (over ocean, “unmixed” pixels only) between the measured and 
model-simulated TOA LW radiances (in Wm 2 sr ] ), between the converted TOA LW fluxes using 
measured LW radiances (both the results from ERBE ADMs and from SDK ADMs are presented) 
and model calculated ones (in Wm' 2 ), both for the nadir-viewing pixels (VZA = 0°) and pixels 
with VZA = 60°. Solid triangles denote the results for clear pixels (data month: May 1994). 

Figure 4b. Same as Fig. 4a, but for July 1994. 

Figure 5a. Same as Fig. 4a, but only for the pixels where model-simulated and observed TOA LW 
radiances agree to with 0.5 Wm' 2 sr'' (data month: May 1994). 

Figure 5b. Same as Fig. 5a, but for July 1994. 
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Figure 6. Scatter plot of TOA LW radiance differences (in Wm' 2 sr'; observation minus 
calculation, nadir-viewing pixels only) against the column precipitable water using (a) TOVS, (b) 
31 atmospheric temperature and humidity profile, along with the corresponding correlation 
coefficient (data month: May 1994). 

Figure 7. Mean and standard deviation of TOA LW radiance differences (in Wm sr’ ; observation 
minus calculation) as a function of ctt c (the standard deviation of cloud top temperatures over the 
16 neighboring ISCCP DX pixels; in K) for the 60° view, ocean ScaRaB pixels classified as high- 
level clouds with optical thickness larger than 1 by ISCCP (data month: May 1994). 

Figure 8. Same as Fig. 4a, but for the “mixed” pixels only. 

Figure 9. Same as Fig. 5a, but for the “mixed” pixels only. 

Figure Al. Zonal mean difference between the ScaRaB DX cloud top pressure retrieved using the 
revised ISCCP cloud program and the program with a coding error in the visible adjustment step 
for May (top panel) and July (bottom panel) 1994. 

Figure A2. Zonal mean difference in the model calculated TOA LW radiative fluxes (in Wm 2 ) 
using May 1994 ScaRaB DX data (revised minus flawed). 
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